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Since several years the preparation and manipulation of a small number of quantum systems 
in a controlled and coherent way is feasible in many experiments. In fact, these experiments are 
nowadays commonly used for quantum simulation and quantum computation. As recently shown, 
such a system can, however, also be utilized to simulate specific behaviors of exponentially larger 
systems. That is, certain quantum computations can be performed by an exponentially smaller 
quantum computer. This compressed quantum computation can be employed to observe for instance 
the quantum phase transition of the ID XY-model using very few qubits. We extend here this notion 
to simulate the behavior of thermal as well as excited states of the ID XY-model. In particular, we 
consider the ID XY-model of a spin chain of n qubits and derive a quantum circuit processing only 
log(n) qubits which simulates the original system. We demonstrate how the behavior of thermal 
as well as any eigenstate of the system can be efficiently simulated in this compressed fashion and 
present a quantum circuit on log(n) qubits to measure the magnetization, the number of kinks, and 
correlations occurring in the thermal as well as any excited state of the original systems. Moreover 
we derive compressed circuits to study time evolutions. 


I. INTRODUCTION 

Simulating a quantum system with a classical com¬ 
puter seems to be an unfeasible task due to the expo¬ 
nential growths of the dimension of the Hilbert space as 
a function of the number of considered systems. This is 
why the classical simulation of quantum behavior is usu¬ 
ally restricted to a few qubits, although the numerical 
methods became very powerful. However, as pointed out 
be Feynman [1] and proven by Lloyd [2] quantum systems 
can be used to simulate the behavior of the other. The 
former being such that constituents can be very precisely 
prepared, manipulated and measured. Many experi¬ 
ments are realizing such a simulation nowadays. Among 
them experiments utilizing ions in ion-traps, NMR or 
atoms in optical lattices (see for instance [3-5] and refer¬ 
ences therein). 

Here we are not concerned about this direct simula¬ 
tion of a quantum system. We are interested in a more 
economical way of simulating certain quantum behaviors. 
To this end, we are using the fact that some classes of 
quantum algorithms, among them those which are based 
on matchgates, can be simulated classically efficiently. 
Moreover, it can be shown that matchgate circuits can 
also be simulated by an exponentially smaller quantum 
computer [6]. There, the classical computation is re¬ 
stricted in space such that the computation has to be 
performed by the quantum computer and cannot be per¬ 
formed by the classical computer. In fact, it has been 
shown that the computational power of matchgate cir¬ 
cuits running on n qubits is equivalent to the one of space- 
bounded quantum computation with space restricted to 
being logarithmic in n [6]. 

As demonstrated in [7, 8] this compressed way of quan¬ 
tum computation can be used to simulated physically in¬ 
teresting behaviours of large systems. To give an exam¬ 
ple, consider an experimental setup, where up to 8 qubits 
can be well controlled. Such a setup can be used to sim¬ 


ulate certain interactions of 2® = 256 qubits. In [7, 8] 
we demonstrated how the adiabatic evolution of the ID 
Ising and more generally the XY-model can be simulated 
via an exponentially smaller quantum system. More pre¬ 
cisely, it is shown there, how the phase transition of such 
a model of a spin chain consisting out of n qubits can 
be observed via a compressed algorithm processing only 
log(n) qubits. The feasibility of such a compressed quan¬ 
tum simulation is due to the fact that the adiabatic evo¬ 
lution and the measurement of the magnetization em¬ 
ployed to observe the phase transition can be described 
by a matchgate circuit. Remarkably, the number of ele¬ 
mentary gates, i.e. the number of single and two-qubit 
gates which are required to implement the compressed 
simulation can be even smaller than required to imple¬ 
ment the original matchgate circuit. This compressed 
algorithm has already been experimentally realized us¬ 
ing NMR quantum computing [9]. In [8] we showed that 
not only the quantum phase transition can be observed 
in this way, but that various other interesting processes, 
such as quantum quenching, where the evolution is nona- 
diabatic, and general time evolutions can be simulated 
with an exponentially smaller system. 

The aim of this paper is to extend the notion of com¬ 
pressed quantum simulation even further. We will con¬ 
sider the XY-model and derive compressed circuits to 
simulate the behavior of the thermal and any excited 
state of the system. To this end, we use the diagonal- 
ization of the XY-Hamiltonian presented in [10]. There, 
the unitary, which diagonalizes the XY-Hamiltonian has 
been derived and decomposed into matchgates. More¬ 
over, it has been shown that this unitary can be uti¬ 
lized to generate the thermal and any excited state of 
the system from a product state. We use these results 
to derive compressed circuit which simulate the corre¬ 
sponding matchgate circuits. They can then be used to 
measure for instance the magnetization of any thermal 
or excited state of the system using exponentially fewer 
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qubits. Apart from that, we also show that the number 
of elementary gates will be smaller in the compressed 
simulation than the direct simulation. 

The outline of the remainder of the paper is the fol¬ 
lowing. First, we will recall the definition of matchgates 
and review how compressed quantum computation can 
be achieved. Moreover, we will recall some basic prop¬ 
erties of the XY-model and its exact diagonalization. In 
Sec. Ill we will recall the results presented in [10]. There, 
the unitary, U, which transforms the Hamiltonian cor¬ 
responding to the XY-model to a Hamiltonian of non¬ 
interacting systems has been derived. As shown there, 
this unitary operator can be decomposed into matchgates 
and can be used to generate the thermal as well as any 
excited state from a product state. In Sec. IV we de¬ 
rive then the compressed circuit corresponding to U. We 
will show that the number of elementary gates required 
for any of those simulations can be even smaller in the 
compressed simulation than the original one. Using this 
circuit it is then easy to derive compressed circuits which 
can be utilized to simulated the behavior of any eigen¬ 
state or thermal state of the model. In Sec. V we will 
then derive compressed circuits for the measurement of 
correlations in the thermal state, which are related to the 
number of kinks in the groundstate, as a function of the 
quenching time and temperature. Moreover, compressed 
circuits which can be employed to measure the magneti¬ 
zation or correlation of the excited or thermal states will 
be presented. Finally we will show how time eolution can 
be simulated with an exponentially smaller system. 


II. PRELIMINARIES 

In this section we first introduce our notation, recall 
the definition of matchgate circuits, and will then review 
how these circuits can be compressed. In the last part of 
this section we recall some basic properties of the XY- 
model. 


A. Notation 

Throughout the paper we denote by X, Y, Z the Pauli 
operators and by 1 the identity operator. The computa¬ 
tional basis will be denoted by \k) for k € {0,1}”. We 
will consider a spin chain of n qubits, where we assume 
that n is a power of 2, i.e. n = 2™ for some integer m. 
For reasons which will become clear later, we will enu¬ 
merate the qubits form 0 to n — I. Moreover, we will use 
the binary notation, j = [jk, •. •, Jo], for j = ji2K 


B. Matchgate circuit and compressed quantum 
simulation 


form A(B B, where A is acting on span {|00), |11)} and 
B on span {|0I), |10)}. Both, A and B are unitary and 
have the same determinant. Denoting by Aij {B^ j) the 
matrix elements of A (B) respectively, we thus have that 
any matchgate is of the form 


GiA,B) 


/All 0 
0 Hi,i 

0 H 2.1 

Va2p 0 


0 Ai2\ 
Bi,2 0 
B2,2 0 ’ 

0 -^2,2/ 


( 1 ) 


with det(A) = det(H). A matchgate circuit of size n is 
defined as a quantum circuit which fulfills the following 
conditions: 


(i) an n-qubit quantum state is initially prepared in a 
computational basis state 

(ii) it evolves under the action of nearest-neighbor 
matchgates, 

(iii) the output is a final measurement on any qubit, say 
qubit k, in the computational basis. That is, the 
output of the computation is the expectation value 

{Zk). 

In [II] (see also [12, 13]) it has been shown that the 
output of any matchgate circuit can be classically effi¬ 
ciently simulated. We will review now how the classical 
simulation can be achieved and will then recall the results 
presented in [6], where an equivalence between matchgate 
circuits and quantum computation running on exponen¬ 
tially less qubit has been shown. 

Let us introduce a set of 2n Hermitian operators 
{a:j}j=o,..., 2 n-i, which satisfy the anticommutation re¬ 
lations 


— 2Sj^k^: (^) 

That is, they define a Clifford Algebra, whose elements 
are any linear combination of products of its 2n gener¬ 
ators Xj. A representation of the generators in term of 
Pauli operators on a system of n qubits is given by the 
Jordan-Wigner (JW) representation defined by 

^2k — Zq . . . Zk—lXk^k+l • • ■ In—1, 

^2fc-t-l — Zq . . . Zk—lYk^k+l • • • Hn —1, 

for fc = 0,..., n — 1. Any matchgate or product of match- 
gates, U, can be written as U = e“*“'^ where H is a 
quadratic Hamiltonian in the operators Xj, that is, 

2n-l 

H = i hj^kXjXk- (4) 

j^k=0 

Here, ft, is a real antisymmetric 2n x 2n matrix [12]. In 
[12] it was proven that for any such unitary, U, it holds 
that 


Let us recall the concept of matchgate and matchgate 
circuits. A matchgate is a two-qubit gate which is of the 


2n-l 

XjU — ^ ^ Bj,kXky 
k=0 


( 5 ) 
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where R = e^“^ G SO{2n) is a real 2n x 2n matrix. 
Using that Zk = —ix 2 kX 2 k+i [see Eq. (3)] and Eq. (5) 
the outcome of the matchgate circuit can be written as 

{Zk) = ( 0 |®” U\-ix2kX2k+i)U\Q)®" 

= Y, R2k,3R2k+U (0|®” {-ix,Xi) |0)®” (g) 

3,^ 

= [RSR^]2k,2k+l, 

where Sj^i = (0|®" {—iXjXi) |0)®" for j ^ I and Sjj = 0. 
That is S' = ]l(8)iU, where 11 denotes the identity operator 
on n = 2™ dimensions. Note that both, R and S are 
2n X 2n matrices. Hence, it is evident from Eq. (6) that 
(Zk) can be classically computed efficiently. 

Let us now briefly outline how such a computation 
can be compressed [6]. Both matrices R and S can 
be interpreted as operators acting on m + 1 qubits, 
which we denote by qubits 0,..., m. Then, S = iY^ 
and the last line in Eq. (6) can be written as {Zk) = 
i (^2k\RYjnR'^\ 2k + l). Simple algebraic manipulation of 
this equation leads to 

(Zfc) =tr[i?^Pi„(fc)i?y^] , (7) 

where pin{k) = \k) {k\ 0 |+y) {+y\m ^ density matrix 

of TO + 1 qubits. The identity in Eq. (7) shows that 
the outcome of the original matchgate circuit, {Zk), can 
be measured by preparing a (to + l)-qubit system in the 
state Pin{k), evolving it under the action of and then 
measuring the last qubit in the Y basis. It is impor¬ 
tant to note here that a decomposition into elementary 
gates of the compressed gate R, which corresponds to the 
matchgates of the original circuit, can be computed with 
a classical computer that is bounded to (!I(logn)-space 
[6]. This ensures that the classical computer behaves 
simply as an encoder and that the quantum computer is 
indeed performing the simulation the original matchgate 
circuit. 


C. XY Hamiltonian 

The ID-XY -Model describes a one-dimensional chain 
of spins, with nearest-neighbor interactions along the X 
and Y directions. Moreover, every spin feels the presence 
of an external magnetic field along the Z direction. The 
Hamiltonian that governs its dynamics is given by 


Cj = {x 2 j +ix 2 j+i) /2. Hence, they can be represented 
in term of Pauli operators as 

Cj = Zo (g) • • • (g) Zj^i (g) (7+ (g) tj + i •••(g) tn-l, (9) 

where = \ {Xj +iYj). Recall that the fermionic an¬ 

nihilation and creation operators satisfy the anticommu¬ 
tation relations {c],c^} = {cj,Ck} = 0, {cj,cl} = Sj^k- 
Clearly, the Hamiltonian H given in Eq. (8) is quadratic 
in the operators c^. It is common to consider the Hamil¬ 
tonian H with JW boundary conditions, which leads to 
the Hamiltonian 


H = H- (10) 

where = {Zq ® ® Zn-i) Xq and Y„ = 

(Zq (g) • • • (g) Z„_i) Yq. We denote by H[c] the Hamilto¬ 
nian H in terms of the c-operators which is given by 

n—1 

Hie] = -Y[9 {cjc] - cjcj) 

3=0 

+ {1 - S) (-CjCj+i + 

+ (1 + <^) + 1 + c]cj+i^ . 

Note that the JW boundary conditions translate into pe¬ 
riodic boundary conditions in terms of the c operators, 
that is c„ = Cq. 

The Hamiltonian H[c] can easily be transformed into 
a diagonal form [10, 14]. This transformation is divided 
into two parts: A Fourier transformation and a Bogoli- 
uvov transformation. The Fourier transformation is de¬ 
fined to map the set of operators {cj}j=o,...,n-i into the 
fermionic operators {fcj}j=o,...,n-i according to 


bk = 


^ n — 1 




3=0 


( 12 ) 


Expressing the c-operators in Eq. (11) by the ^-operators, 
the Hamiltonian takes the form 

n— 1 

H[b] = -Y + i/3j {b-jb^ + bl^bf^ + g . 

3=0 

(13) 

Here and in the following we used the definition b-j = 
bn-j, for j = 0,..., n — 1. The coefficients aj and l3j are 
given by 


H = - 


n—1 n—2 

'Y 9^3 + Y^ i^3^3 + i + ^XjYj+i) 

3=0 j=0 


( 8 ) 


The parameter g represents the ratio between the 
strengths of the external magnetic field and the near¬ 
est neighbor interactions, while the factor 6 € [0,1] rep¬ 
resents an anisotropy between the X and Y direction 
of the spin-spin interactions. The Jordan-Wigner trans¬ 
formation [14] can be used to transform the spin oper¬ 
ators into the fermionic operators, {cj}j=o,...,n-i, with 



(1 -I- (5) — g, and 


/3j = sin 



{ 1 - 6 ). 


(14) 


The diagonalization of the Hamiltonian is completed 
by a Bogoliuvov transformation that maps the oper¬ 
ators {bj}j=o,...,n-i into a set of fermionic operators 
{aj}j=o,...,n-i which are defined as 

Gj = Ujbj - ivjblj 


(15) 
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Here, the coefficients Uj = cos j and Vj 
are functions of the angle 



9j = tan ^ 



(16) 


Expressing the Hamiltonian in terms of the a-operators 
one finds 


n—1 

H[a] = Y^ ej 

The energies are determined by 

= -2 [aj cos {ej) + Pj sin (%)], (18) 

e.g. the groundstate energy is Eq = 

Eq. (16) defines the angles 9j up to the addition of tt, the 
sign of the energies ej is not defined uniquely. In order to 
have a consistent labeling of the energies for every value 
of g, that is, e.g. that eo is the gap between the first 
exited state and the groundstate energy for any value 
of g, one imposes that the energies tj have to be non¬ 
negative for all j. This constraint determines uniquely 
the value of 9j. Note that this choice amounts to keeping 
the same ordering of the energies, i.e. the first excited 
state corresponds to the energy Eg -I- cq. 


- 2 


(17) 


1. Phase transitions 

For any constant 5 the XY Hamiltonian given in 
Eq. (8) exhibits a quantum phase transition. One distin¬ 
guished two different scenarios at the phase transition: 
either the spectrum exhibits a level crossing for a finite 
system size, which occurs if the energy of the ground- 
state is degenerate at the critical point; or the ground- 
state energy is not degenerate for any finite system size, 
but the energy gap between the groundstate and first- 
excited state tends to zero in the thermodynamic limit. 
This situation is referred to as an avoided level crossing. 
In both cases the groundstate energy is non-analytic at 
the phase transition, which is reflected as an abrupt or 
discontinuous behavior of macroscopic observables, such 
as the magnetization, at this point. The relevance of 
quantum phase transitions relies on the fact that they 
occur at zero temperature and are caused solely by the 
quantum-mechanical nature of the system under consid¬ 
eration. 

The case we consider in this work is the XY Hamil¬ 
tonian with Jordan-Wigner boundary conditions. Under 
these boundary conditions, the system exhibits a level 
crossing at the critical point gc- This is due to the sym¬ 
metries of the Hamiltonian with respect to parity and 
momentum [8, 15]. The spectrum of any subspace with 
fixed parity and momentum presents no phase level cross¬ 
ing. However, in any of these subspaces a level crossing is 
still present when considering the thermodynamic limit. 


In Sec. V A 1 we present a compressed quantum circuit 
that can be used to measure the magnetization of the 
groundstate of the Hamiltonian displaying a phase tran¬ 
sition (see also [8]). This circuit is extended in Sec. V 
where we present a compressed quantum circuit that can 
be used to measure the magnetization of any thermal and 
excited state of the Hamiltonian. 


2. Quantum Quenching 

A relevant physical evolution that can be investigated 
in systems with quantum phase transitions is denomi¬ 
nated quantum quenching. That is, to observe the evo¬ 
lution a the system prepared in the groundstate of the 
Hamiltonian for a certain value of g when it is driven 
through the critical point by changing the parameter g 
abruptly. In the cases where the Hamiltonian does not 
exhibit a level crossing and the driving speed is slow 
enough, the system remains in the instantaneous ground- 
state of the system (adiabatic evolution). However, when 
the quenching time t is too short, the evolution gener¬ 
ates excitations in the system, these excitations can be 
quantified as a certain amount of defects v in the qubit 
alignment, commonly referred to as kinks. The number 
of kinks, therefore, depends on the quenching speed, and 
tends to zero when t —>■ oo. The form of the dependence 
of the number of kinks on the quenching speed is pre¬ 
dicted by the so called Kibble-Zureck mechanism [16, 17], 
that predicts a dependence of the form v oc . The 
same dependence has been observed for systems of differ¬ 
ent nature that also present a quantum phase transition 
(see for instance [18] and references therein). 

In Sec. V A 2 we present a compressed quantum cir¬ 
cuit that can be used to measure the number of kinks 
in the Ising Hamiltonian [(5 = 0 in Eq. (8)] as a func¬ 
tion of the quenching speed. There we also compute the 
dependence of the correlations, which are related to the 
number of kinks, on the quenching speed for the case of 
finite temperature. 


III. MATCHGATE CIRCUIT 

In the previous section we have seen that the diago- 
nalization of the Hamiltonian is achieved by transform¬ 
ing the set of fermionic operators {ci}i=o,...,n-i into a 
different set of fermionic operators {ai}i=o,...,n-i- Let 
us denote by , the corresponding transformation on 
the fermionic operators. In this section we will review 
the construction of the corresponding unitary matrix U, 
that transforms the Hamiltonian as 

UH[a]U^ = H[c]. (19) 

Then, we will recall now the decomposition of U into 
matchgates presented in [10] [20]. 

First, note that due to the JW transformation, one can 
associate to an n-qubit state I'k) a fermionic state in the 
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following way, 

IV’[a])= V'Nio.-.i,.-! l*0,---,*n-l) 

= ( 20 ) 

Here, and in the following |H[a]) represents the vacuum 
state of modes at, that is, at |H[a]) = 0, Vt. Note that us¬ 
ing this notation the operator aj is associated with qubit 
j. As the fermionic operators do not commute, any re¬ 
ordering of the operators might result in a change of signs 
in Eq. (20). Sometimes, however, it is more convenient to 
work with a different order, especially as we are seeking 
for a decomposition of U into nearest neighbor match- 
gates. We will explain in the following how to take the 
phases into account. 

In order to construct a nearest neighbor matchgate 
circuit U, the transformation U is decomposed into el¬ 
ementary transformations, that act non-trivially only on 
2 fermionic operators, as we will describe in detail in the 
following subsections. These transformations are con¬ 
jugated by reordering transformations that permute the 
position of the qubits. Due to the anticommuting rela¬ 
tions of the fermionic operators, swapping, e.g. qubit j 
and k ^ j, yields a phase (—K*fc-i) jjj 
each of the terms of Eq. (20). It is convenient to intro¬ 
duce a vector /x = (/xq, • ■ ./i„_i), that is a permutation 
of the vector (0,..., n — 1), whose components indicate 
which fermionic operator is associated to which qubit. 
This is done by generalizing the notation of Eq. (20) to 

( 21 ) 

where denotes here an arbitrary set of 

fermionic operators. Let us denote by Sj^k the fermionic 
swap operator, which maps aj to and vice versa. 
The action of the transformation Sj^k corresponds to the 
permutation of the two components of /x, /ij and /Xfe. 
Stated differently, we have iply, coincides with 

(0,..., n — up to the sign, which re¬ 

sults from the permutation. 

In the following section we will consider transforma¬ 
tions that act non-trivially only on two fermionic opera¬ 
tors, associated to consecutive qubits, e.g. 

In those cases the corresponding gate acts nontrivially 
only on qubits 21 and 21-1-1 [21]. The correspondence 
between the computational basis states of qubits 21 and 
21-1-1 and the fermionic operators is then [see Eq. (21)], 


\[ki,k2][y,fi])2i ,2/+l 


{yLf' {vL+i) I^M) 


21,21^1 ’ 


( 22 ) 


where, the vacuum state |D[t/]) 2 ; 21+1 represents the vac¬ 
uum state of modes and y^j, 2 i+i ■ 

In the previous section we have followed the conven¬ 
tional procedure used in the literature to diagonalize the 
XY Hamiltonian. There, the diagonalization has been 
achieved by applying a Bogoliuvov transformation, B, 
and a Fourier transformation, T. In the following subsec¬ 
tions, we recall how to construct the matchgate circuits 
corresponding to each of these two transformations [10]. 


A. Bogoliuvov Transformation 


The Bogoliuvov transformation B maps the fermionic 
operators {ai}i=o,...,n-i into the operators {&i}i=o,...,n-ij 
according to Eq. (15). It is easy to see that B can be writ¬ 
ten as B = ri/Lo where Bi transforms the operators 
{ai,a-i} into the operators {bi,b-i} [see Eq. (15)], for 
I = 1,..., n/2 — 1 and Bq maps {oq, a^} into {&o, }■ 

As the two modes which are mixed by the transforma¬ 
tion Bi, have non consecutive indexes we need to reorder 
the qubits, by applying a permutation iSeog to obtain 
a product of gates acting on consecutive indices. This 
transformation permutes the position of the qubits in 
such a way that the qubits I and —I = n — I, and 0 
and ^ become consecutive before the transformation Bi 
acts, and are returned to the original position afterwards. 
As was explained above, we can define a vector /x that 
indicates the index of each fermionic operator after the 
reordering transformation [see Eq. (21)]. We have chosen 
a particular transformation which corresponds to 

the mapping of the vector (0,..., n — 1) to the vector 


A^°s= (0,|,2,-2,...,2Z,-2;,...,^-2,-|+2, 

Ti Ti \ 

- - 1 , -- + 1 ,..., 2Z - 1 , -2Z + 1 ,..., 1 , -Ij. (23) 

Using the following notation for 0 < ? < f — 1 


ri 


I for I even, 

for I odd. 


(24) 


the components of satisfy that (Af°®, = 

(1,-1) for 1 < / < f-1 and (Ao°®, A®°®) = (0, f). There¬ 
fore, it is clear that the transformation Bi now acts non- 
trivially on the pair of operators \ a,Bog, a,Bog L which 

L ^2ri ^2ri+l ) 

are associated to the consecutive qubits 2r; and 2r/ + 1 
[see Eq. (21)]. Note that both indexes I and n are integer 
numbers in the interval [O, f ~ l]; the index ri indicates 
the position of the pair (I, —1) in the vector A^°®. 

Let us now compute the gate Bi, which corresponds 
to Bi- First note that the transformation Bi, maps 
the vacuum state | 0 [a]) 2 r, 2 ri+i ^ vacuum 

state |H[5])2 ,,j 2 ri+i- order to compute this new vac¬ 
uum state, note that ai |D[a]) 2 ,ri 2 ri+i ~ ^ ^ ^ 
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A 2 ^^i}- Hence, one can write the projector onto 
the vacuum state as 




2 n, 2 n+i “ 

2ri 


(25) 


Expressing the operators a,Bog and a,Bog in the b op- 

''2r, ''2r,+l 

erators [see Eq. (15)], the vacuum state | 12 [a]) 2 r, 2 ri+i i® 
mapped to the state 


+ ivib 


X Bog 


|11[5]) 


2r,,2r, + l 


(26) 


It is now straight forward to compute the two-qubit 
gates Bi that correspond to Bi by noting that Bi maps 


[ki,k 2 ] a,A®°® ^ to [ki,k 2 ] 6 , A®°® ^ for ki,k 2 


{0,1}. The relation between the elements of these two 
bases can be computed using Eqs. (22), (26) and (15), 
leading to 


[fci,fc2] [a,AB°s]\ 

/ 2ri,2ri + l 


= ( + ivib. 


-'2^ 


2r-i+l 


Ulbleag - iVib,Bae 


Uil + iviblBagbK^g 


>DOg“X 

'^2ri+l 


|H[&]) 


2ri,2ri + l ' 


(27) 


By evaluating Eq. (27) for k 2 ,ki € {0,1} and us¬ 
ing Eq. (22) on the right-hand side of Eq. (27) to 


identify the basis states 


[j1,J2][6,A«°8] 


one 


2ri,2ri + l 

can compute the components of jj] [fc^ fc 2 ] “ 

([ji,j 2 ][&, A®°®] [fci,fc 2 ][a, A®°®]y Doing so, one can eas¬ 
ily find that 


Bi = 


/m; 0 0 ivA 
0 10 0 
0 0 10 
\ivi 0 0 uiJ 


(28) 


As uf +vf = 1, this matrix is a nearest neighbor match- 
gate that acts non trivially solely on qubits 2 r; and 2ri+l. 
Recall that ui and vi depends on 0;, which has to be cho¬ 
sen such that e; > 0 [see Eq. (18)]. That is, depending on 
g the the angle is either 0; as given in Eq. (16) or 9i + tt. 

The matrix Bq can be computed similarly, but it 
requires special attention because as can be seeen in 
Eq. (15), oo = uobo - ivobl and = uiibii - iviib\. 
That is, the the transformation Bq mixes the operators 
with same indexes, which clearly leads to a different gate. 
In order to compute it let us first note that the coef¬ 
ficients aj and (3j dehned in Eq. (14) take the values 
ao = ^ + 5 — g and /3o = 0 for j = 0. Therefore, the angle 
6*0 = tan“^ (^) values 0 or tt. Using that 

the energy eg is given by cq = — 2(1 -|- (5 — 5 ) cos( 6 * 0 ) and 


that 6*0 is chosen such that eo > 0 , we have 


70 = 


0 for g > 1 -I- (5, 

TT for g < 1 -I- (5. 


(29) 


It follows from this expression and from Eq. (15) that 
(ao,af) = ibo,b^) ior g > l-h(5 and ( 00 , 0 ^) = -i(&t, 6 j) 
for g < 1 -|- (5. Thus one can easily compute that 


Bo = 


1 

Xq 0 Zi 


for g > 1 + S, 
for g < 1 + S. 


(30) 


Note that if g < 1 -I- 5 the gate Bg is not a matchgate as 
the condition on the determinants is not satisfied. How¬ 
ever, in Sec. IVA, when compressing the gates Bi, we 
will show that even in that case a corresponding com¬ 
pressed gate can be found. Note that the appearance 
of the operator Xg is due to the level crossing of eigen¬ 
states with even parity with eigenstates with odd parity. 
For instance, the groundstate has even parity in the re¬ 
gion where g > 1 -\- 5 and odd parity in the region where 
g<l + 5 [ 22 ]. 

The circuit B corresponding to the transformation B 
is now obtained by multiplying the gates Bg given in 
Eq. (30) and Bi in Eq. (28), that is 


B = Bro0---0Br^_^. (31) 


B. Fourier Transformation 

The Fourier transformation X maps the operators 
{fej}j=o,...,ra-i into the operators {cj}j=o,...,n-i, where 


Cfc = 


n—1 

E 


exp 


—i27Tjk 


Uj. 


(32) 


Below, we review the procedure presented in [10] to write 
the transformation as a concatenation of transforma¬ 
tions which involve only two fermionic operators at a 
time. As each of those elementary transformations cor¬ 
responds to a matchgate, we will obtain in this way the 
matchgate circuit corresponding to B. 

The idea behind the decomposition of B is that the 
mapping of the operators {fcj}j=o,...,n-i into the oper¬ 
ators {cj}j=o,...,n-i can be accomplished by applying a 
sequence of maps B^^\ for 0 < s < m— 1 , to the operators 
bj, i.e. B = J7(in-i) o • • • o J 7 (o)_ transforma¬ 

tion, J-(o), maps the fermionic operators {6j}j=o,...,n-i 
into a different set of fermionic operators {a;j^^}j=o,...,n-i, 
which is defined in such a way that it always mixes only 
the pair of operators {bj, bji}, where the indices j and j' 
differ only in the {m — 1 )*^ bit of their binary represen¬ 
tations. This procedure is repeated iteratively, in such a 
way that at the s*'' step, the indices of the pair of oper¬ 
ators that are mixed differ only in the {m — s — 1 )*^ bit 
value. Below we explain the procedure in detail. 
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In the following we will refer to j = either 

as the number or its binary representation [jm-i, ■ ■ ■, Jo] 
whenever it does not lead to any confusion. The phase 
in Eq. (32) can now be expanded as 




n 


a-l-l 


(33) 


1}, we find 


Ck 


E 

io=o 


1 

Vi 


e 








1 


(s 

[ko,...,ks-i 


,3m 


-i,---,3o]' 


(38) 


Let us denote the initial operators bj = ^ ,Vj, 

to use a notation compatible with the iterative process 
that we are going to utilize. Using the decomposition of 
the phases given in Eq. (33), Eq. (32) can be written as 






inkjQ 


io=o 


e 2 ’ 




-i7rfcj^_i (0) 

^Ur.-U...,j0V 

(34) 

Computing only the sum over the index jm-i we obtain 


J„_i=0 


V — 1 

^ V2 


jm-l=0 


bm-l,..-Jo] 


V2 L' 
- ..(1) 


(0) 

'^[Oj™_2....Jo] 


i_li'=o.y.(0) 

+ 1 C 


= X 


[fco Jto- 2 ....Jo]’ 


(35) 


where fco is the last bit value of k. Recalling that 
[fco) Jm-i, • • •) Jo] corresponds to the number fco2'"“^ + 
one can see that Eq. (35) defines a set of n 
fermionic operators {a;(^^}i=o,...,n-i as a function of the 
operators {a;(°^}i=o,...,n-i- Explicitly, we have that 


(1) 

lOjm-l,. 

■■Jo] 

I 

“Vi 

_^[0j™-2,^ 

+ 

■■Jo] ]1 Jm- 2 ,^^^ Jo]_ 

(1) 


1 



]1 Jm-l,. 

■■Jo] 

“Vi 

_^[0j™-2,^ 

_ 

■■Jo] ]1 Jm- 2 ,^^^ Jo]_ 


Let us call the transformation that maps the original 
fermionic operators {a;^°^}i=o,...,ra-i into the fermionic 
operators {x-^^}i=o,...,n-i) according to Eq. (36), 

Note that the transformation mixes pairs of opera¬ 
tors whose indices differ only in the {m — 1)*^*' bit value. 
Replacing Eq. (35) in Eq. (34) we obtain 




E 

io=o 



2^ — 1 




jm.-2=0 


™(1) 

■^[fcoJm- 2 ,... Jo]' 


(37) 


Similarly as before, we define the operators 
{a;(^~''^^}i=o,...,n-i as the sum over the index jm-s-i in 
Eq. (38), that is 


X 


(s+l) 

]fc 0 ,...,fcs Jm-s- 2 ,... Jo] 


V2 




,C) 


(- 1 ) 


ks ^-i7T(p{k,s) ^ 


[ko,...,ks-l,0,jm-s-2,---,jo]~ 

,»])■ P9) 


]fc0,...,fcs-l,l Jm-s-2, 


In this general case, the phase difference be¬ 
tween the operators k i o j 2 Jo] 

fe 1 1 j 2 Jo] additional component 

(p{k,s). It can be extracted by evaluating explicitly 
e 2 = for jm-8-i = 1, that is 


e-^ = 2 ' ^ (40) 

with (p{k, s) = X]/*=o In the last equality we used 

the fact that the terms in the sum with I > s lead to 
trivial phases. In the same way as was defined via 
Eq. (36), we define the transformation as the one 
that maps the fermionic operators {x(*^}i=o....,n-i into 
the new fermionic operators given in 

Eq. (39). The transformation mixes pairs of operators 
whose indices differ only in the (to — s — I)*^ bit value. 

The last transformation of the succession 
••• o J-C), i.e. maps the fermionic oper¬ 
ators ^^}i=o,...,n-i into the fermionic operators 

Til® output of this transformation coin¬ 
cides, up to a permutation of the indices, with the desired 
set of operators {ci}i=o,...,n-i- In fact, setting s = to — I 
in Eq. (38) we obtain 

C[fc„_i,...,fco] ~ ^lko]...,k^-i]' (^V 

Let us now have a closer look at the transforma¬ 
tion for an arbitrary s. Due to Eq. (39) we have 
that the transformation can be written as = 
o- ■ -o where each map for 0 < / < | — 1, 

acts non-trivially only on one pair of operators, namely, 
|j where r; ^ corresponds to the num¬ 
ber whose binary representation is 


[^m—25 ■ ■ ■ 5 ^m—s—ly O7 ^m—s — 2 i ■ • ■ : ^o] ■ 


(42) 


Repeating this procedure s times, for s € {0,..., to — 


Hence, for each I, the indices of the two operators that 
are non-trivially transformed, differ by 2™“^“^, i.e. by 



their [m—s—lY^ bit value in their binary representation. 

(s) 

More explicitly we have that J~i maps the pairs of op¬ 
erators fhe pairs of operators 


{= 


s+l) ^(«+l) 

<.» >‘^r,.+ 2 " 


_i I according to 




„ — rr \ , ~ ^ ^_|_ 2 '^-s-l 


X 


(«+l) 

n,s+2^ 




I ^(s) _ p-™v(l,s) ^(s) 


where the phase is given by 

s-l 


j=0 


(43) 

(44) 


Note that this is just a redefinition of the phase (f{k, s), 
given in Eq. (40), that takes into account that the binary 
bits of I in Eq. (42) are sorted differently than those of k 
in Eq. (38). 

Now we are in the position to review how to construct 

the matchgate associated to the elementary transforma- 
(s) 

tions Ti , yi, s. Let us assume that before and after each 
of these operations, appropriate reordering operations are 
implemented such that the ordering of the fermionic op¬ 
erators is described by a vector [see Eq. (20)], with 

, for 


components and = r/^g -I- 2 

I = 0,..., I — 1. Therefore the transformation maps 

now the operators < f into 

( ^ 2 i ^21 + 1 j 

that are associated to consecutive qubits 21 and 21-1-1. 

(- 5 ) 

Hence, the gate corresponding to Jq , which we de¬ 
note by acts non trivially only on the 21*"^ and 

{21 + 1)*^ qubits. Unlike the Bogoliuvov transforma¬ 
tion, this transformation is passive and therefore does 
not modify the vacuum state, i.e. = |H[a;i*i]) 

[23]. The two-qubit gate Fi that corresponds to the trans¬ 
formation Fi can be computed by noting that Fi maps 

■[fci,fe][x(*),A(*)]) to [ji,J 2 ][x(^+'),a(^)]) 

/ 21 , 21+1 / 21,21 + 1 

Using the notation given in Eq. (22) with operators 

and Eq. (43) one obtains that 


(s) 


')m —s— 1 
(^) 


X (3) 




[fcl,fc2]M^\AW] 


21,21+1 


(s+l)t I (s+l)t 


X{3) 


x(-) 


^/2 


ki 


T(p{l,s) 


■ ' ' (,) I 

■^21 + 1 / 


X(3) 


72 


k2 




21,21+1 


(45) 


By evaluating this expression for ki,k 2 G {Oj l}j 
and using Eq. (22) to identify the basis ele¬ 


ments 

hand 

pute 


on the right- 

side of the equation, one can com- 
the matrix elements [^i][ii 02 ].[fci.fe 2 ] 


([ji,J2][7^+i\a(^)] [fcl,fc2][7^^A(^)]). 


One can 


(s) . . 

easily verify that, F^ , which is a matchgate acting non 
trivially only on qubits 21 and 21-1-1, is given by 


/I 

0 

0 

VO 


0 


cx{l,s) 

' w 

a{l,s) 


V2 

0 


0 

1 

V2 

1 

V2 

0 


0 \ 
0 

0 

-a{l,s)J 


(46) 


where a{l,s) = e^'^vi.ks) ip{l,s) given in Eq. (44). 
As is decomposed as a sequence of transformations 
Fi ' ioT 0 < I < f — 1, the matchgate product associated 
to is then given by 


=Fo^®^(g)---(g)Fl®y. (47) 

The matchgate circuit corresponding to the whole Fourier 
transformation is then given by the product of circuits 
F^^'>, conjugated by appropriate reordering matchgates 
that correspond to the reordering transformations dis¬ 
cussed above. 


C. Reordering operations 

As mentioned before, the different steps of the trans¬ 
formation which diagonalizes the Hamiltonian have to 
be conjugated by reordering transformations. The rea¬ 
son for that is that in order to obtain nearest-neighbor 
matchgate circuits corresponding to the Bogoliuvov and 
Fourier transformations, only fermionic operators in con¬ 
secutive positions have to be mixed at any step of the 
transformations. The transformation U is thus given by 
[24] 

U = 5 (™) o ... o 

(48) 

where the transformations and for 0 < fc < m 
represent the reordering transformations. Obviously, any 
reordering transformation can be accomplished by a se¬ 
quence of reordering transformations that corre¬ 

spond to the swapping of qubits I and I + 1 respectively. 
In Appendix A we derive the explicit form for each of the 
permutations. Here, we just recall that the gate, 
corresponding to which can be computed as be¬ 

fore, corresponds to the fermionic swap, i.e. 

/I 0 0 0 \ 

Q _ 0 0 1 0 

^1,1+1 - 0 1 0 0 
Vo 0 0 -1/ 


( 49 ) 
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The gate Sj^k, corresponding to Sj^k which exchanges 
with is then given by 

fc-i 

Sj,k = n (50) 

i=j 

The matchgate circuit U can now be written as 


We consider hrst the case g > 1 + 6. In this regime 
Bo = 1 [see Eq. (30)] and therefore we can extend the 
definition of Bi in Eq. (28) to the case 1 = 0. As explained 
in Sec. Ill A, the gate Bi acts on qubits 2r; and 2r/ + 1, 
where r; is given in Eq. (24). It is easy to see that the 
gate Bi can be written as 

Bi = (53) 


JJ = g(m) p(m-l) (51) 

where is defined in Eq. (47), B in Eq. (31) and the 
gates S are given in the Appendix A as a function of the 
fermionic swaps Sj^k- In the next section, we derive the 
compressed gate corresponding to U using Eq. (5). 


IV. COMPRESSED CIRCUIT 

As shown in Sec. Ill, the matchgate circuit U is de- 
composed into matchgates B;, ’ and for I G 

{0,..., I — 1} and j,k G {0,..., n— 1}. In this section we 
construct the compressed gates corresponding to each of 
these gates, namely Rbh Rfi and Rs^ This is done by 
following the procedure described in Sec. II. That is, we 
write an arbitrary matchgate G in the form G = e“*^°, 
with Hq of the form of Eq. (4), from where we can ex¬ 
tract the matrix of coefficients ha- The corresponding 
compressed gate is then given by Rq = . 

Although Rb,, RpM and Bs-. are all the compressed 
gates required to construct the compressed circuit R, it 
is possible to use much less gates by taking the symme¬ 
tries of the Hamiltonian into account. As we will see, 
the compressed gate corresponding to the whole prod¬ 
uct of matchgates given in Eq. (47) 

can be more economically implemented than the prod- 
uct of the compressed gates corresponding to , Vh 
The same applies to any of the factors of the product 
given in Eq. (51). For this reason, we write down the 
compressed gate corresponding to U, as the product 

R = Rgiin) Rp(m-l) Rg(rr,-1) . . . Rp(0) Rg (0) RgBag Rb Rs^°e 1 

(52) 

where we have used the notation Ry to denote the com¬ 
pressed gate corresponding to a unitary V. We omitted 
in Eq. (52) the dependence oi R on g and 6. However, 
whenever we want to stress this dependence, we will write 
R{g,S) for R given in Eq. (52). Below we are going to 
provide the expression for each of the factors. 


A. Compression of the Bogoliuvov transformation 

In this subsection we derive the compressed gate, 
Rb, corresponding to the Bogoliuvov transformation B, 
which is a product of the matchgates Bi [see Eq. (31)]. 
In order to do so we compute Rbi, ke. the compressed 
gate corresponding to Bp 


where 9i is given in Eq. (16). Using the Jordan-Wigner 
representation given in Eq. (3), Eq. (53) becomes B/ = 
, with a quadratic Hamiltonian 


Hbi 



{x4ri^4ri+3 F X4ri-\-lX4ri+2 


h.c.), 


(54) 


where h.c. denotes here and in the follow¬ 
ing the hermitian conjugation of the preceding 
terms. Writing iJs, = * ajjXfc, with 

hB, = -f- (|4n) (4r/ -1-31-1- I4n -I- 1) (4n -I- 2] - h.c.) one 
can easily compute the compressed gate , 

which is then given by 



Here, Bs, acts non-trivially only on the basis elements 
{I4r/),..., I4r/-I-3)}. Note that Bs, = 
where Bb, is diagonal and the structure of Vbi is in¬ 
dependent of I (see below). Due to the independency of 
Vb, on I (only the basis elements on which it acts non- 
trivially depend on 1), Rb can be written as 


Bb = n = vIDbVb. (56) 

Here, the gate Bb is diagonal and is given by 


^B=W Dbi 
1=0 

= E |0) (0|+e*^Ml)(l|) 

1=0 

and the gate 

„_i /z 0 0 U 

y, VT TA 11 ijOzlO 

ub= n =1®^ 0 -i 10 

\-z 0 0 1/ 


(57) 


(58) 


acts non-trivially only on the last two qubits and it is 
independent of the angles Op 

In the regime where g < 1 -I- d Bg is given by Ag 0 Zi 
[see Eq. (30)]. Even though it cannot be written as the 



10 


exponent of a quadratic Hamiltonian one can neverthe¬ 
less find an expression similar to Eq. (5) as we will show 
in the following. Note that Bq commutes with Xj for 
j = 0,2,3 and anticommutes with any other Xj. Thus 
we have that 


2n-l 

k=0 


(59) 


with 


i?B„=2(|0)(0| + |2)(2| + |3)(3|)-l. (60) 

In contrast to the matrix R in Eq. (5), the matrix Rbq 
is orthogonal but its determinant is not one. However, 
Eq. (59) ensures classical simulability and compressibility 
of the matchgate circuit, independently of the determi¬ 
nant of Rbo [25]. 

In summary, for g > I -|- <5, the compressed gate cor¬ 
responding to the Bogoliuvov transformation is given by 
Rb given in Eq. (56), while for (/ < 1 -|- d, it is given by 
the product RbRbq- 


B. Compression of the Fourier transformation 

In this section we compute the compressed gate R„(s) 

(s) 

corresponding to the matchgate F; . The compressed 
circuit corresponding to is then given by the product 
= n;=o Rp(s). Below we derive Rp(s) which is 
computed similarly to Rb in Eq. (56). 

The gate given in Eq. (46) acts non-trivially on 
qubits 21 and 2/ -|-1. The corresponding compressed gate 
can be easily found to be 


is diagonal and the gates 


Fg = I]^ Rgi ^ -I 


1=0 



vf = n r. 


(») 


1=0 


/I 0 0 0 \ 

'o 1 0 o' 
n 0 — — 

^ ^ V2 V2 


and 


(64) 


act non-trivially only on the last two qubits. 


C. Compression of the reordering steps 

As mentioned in the previous section, every reorder¬ 
ing step consists of a product of permutations that are 
associated to a product of the matchgates Sj^k given in 
Eq. (50). Below we provide its corresponding compressed 
gate Rsj,k- The compressed gates in Eq. (52) that cor¬ 
respond to each reordering step can all be computed as 
a product of matrices RSj^k detailed circuits see 

Appendix B). 

The matrix Sjj+i [see Eq. (49)] can be written, up to 
a global phase, as 

C. . — „p(0C2jX2j + l+X2j + 2X2j+3+X2j + lX2j + 2-0C2jX2j + 3-'tL-C.) 

(65) 

where the operators Xj are given in Eq. (3). Hence, the 
corresponding compressed gate is given by 

Fs,„+i = (l -b 1 j) {j -b 1] -b \j + 1) {j\ 

^ ( 66 ) 

- |j) (il - \j + 1) (j + 11)® tm, 


R 


,(») 



(61) 


where F„(s) is diagonal and the structure of the matrices 
Rgi and Vp(s) does neither depend on I nor on s. The 
matrix R„(s), and therefore all the elements in the de- 
composition given in Eq. (61) act non trivially only on 
the basis elements {[4?) ,\4l + 3)} [see Eq. (5)]. Com¬ 
bining all that leads to 


RpG, = Yl RpM = RgVpDpg)VI (62) 


where the operator in parenthesis is acting on the qubits 
0,..., m — 1. As Sj^k can be easily obtained by multiply¬ 
ing consecutive fermionic swap gates, the com¬ 

pressed gate, Rs^ can be obtained by multiplying con¬ 
secutive gates of the form of Eq. (66). Explicitly we have 

{k\ + \k) {j\-\j) {j\-\k) {k\)(^lm. (67) 

In the subsequent section we determine the number of 
elementary gates required to implement the compressed 
circuit. 


D. Comparison of the number of elementary gates 


Here, the gate 

Dpis, = \0) (0l^_i®l+ll)(ll^_i® 

E [l0(/|®(e-*"‘^('’^'|0)(0L+F-^('’^)ll)(ll)' , 

1=0 

(63) 


As mentioned at the beginning of this section, the sym¬ 
metries of the Hamiltonian allow to decrease the num¬ 
ber of single- and two-qubit gates required to implement 
the compressed gate F compared to the number of near¬ 
est neighbor matchgates required to implement U in the 
original matchgate circuit. While the decomposition of 
U into matchgates considered here requires gates 
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[10], the number of single- and two-qubit gates that suf¬ 
fice to implement R scales as 0{nm) = 0[n\og{n)]. In 
order to see this, we determine the number of elementary 
gates of each factor of the decomposition of R given in 
Eq. (52). 

As we will show in Appendix B, for each 1 < s < m — 1 
the matrix Rg(s} is a single two-qubit gate. Given that 
i? 5 (o) can be decomposed as a product of m — 1 Rs(s) 
matrices [see Eq. (119) in Appendix B], it follows that 
Rg(o) can be implemented with m — 1 two-qubit gates. 
The gate Rgim) can also be decomposed into a similar 
product of [to/ 2J two-qubit gates [see Eq. (121) in Ap¬ 
pendix Bj. Hence, the implementation of Rgm, for any 
k, requires at most 0{m) two-qubit gates. The matrix 
i? 5 Bog can be implemented with 0{w?) = 0[\og{n)‘^] ele¬ 
mentary gates as can be seen as follows. Due to Eq. (124) 
(see Appendix B) it is clear that Rgsog can be written 
as a product of three gates controlled by two qubits plus 
some small constant number of elementary gates. From 
the decomposition given in Appendix B it can be seen 
that each of these gates require at most 0{m?) elemen¬ 
tary gates. 

The implementation of Rp(s) can be done with 0{s) 
single- and two-qubits gates for 0 < s < m—1. In order to 
see that, first note that the matrices Rq and Vp appear¬ 
ing in the decomposition of Rp(s) given in Eq. (62) are 
single two-qubit gates. Additionally, using the expres¬ 
sion of ip{l,s) given in Eq. (44), the gate Dp(s), given in 
Eq. (63), can be written as a product of 2s single-qubit 
gates controlled by two qubits. Therefore, Dp(s) can be 
implemented with 0{s) elementary gates [19]. It follows 
that 0{rn?) elementary gates are sufficient to implement 
all the compressed gates associated to the Fourier trans¬ 
formation. 

The matrix Rp can be implemented with 0{nm) 
single- and two-qubit gates, which is the leading scaling 
order in the number of elementary gates that decompose 
R. This scaling is due to the implementation ot Dp [see 
Eq. (57)]. One can easily see that Dp can be decomposed 
into a product of | — 1 single-qubit gates controlled by 
m — 1 qubits, each of which can be implemented with 
0{m) elementary gates [19]. Note that due to the in¬ 
direct dependance of the phases 9i with respect to the 
index /, it seems unfeasible that Dp can be implemented 
in a more efficient way. Hence, the compressed circuit 
can be implemented with 0{mn) elementary gates. 

Let us point out that the implementation of the 
gates corresponding only to the Fourier transfor¬ 
mation, i.e. the implementation of the product 
Rg{m.)Rp(m-i)Rg(r„.-i) ... Rp(o)Rs(o) cuu bc done with 
O [log(n)^] single- and two-qubit gates. In contrast to 
that, the number of nearest neighbor matchgates re¬ 
quired to implement the Fourier transform in the original 
matchgate circuit scales as 0{n'^) [10]. 

In the subsequent section we are going to utilize the 
compressed gate R, given in Eq. (52), to construct com¬ 
pressed circuits which can be employed to measure phys¬ 
ically relevant quantities using exponentially less qubits. 


V. APPLICATIONS 

In Sec. HI we recalled how to decompose the matrix U, 
which is diagonalizing the XY-Hamiltonian, into a prod¬ 
uct of matchgates. Any excited state or thermal state 
can be generated by applying U to a product state, as 
we will recall in this section. Once this state is prepared, 
it is possible to measure observables whose expectation 
values reveal physically interesting properties of the sys¬ 
tem. The matchgate circuits consists of three steps: the 
preparation of a factorizable initial state, the evolution 
of the system by the unitary U, and the measurement 
of the desired observable. The aim of this section is 
to provide compressed quantum circuits that simulate 
these matchgates circuits. The condition that the given 
matchgate circuits are compressible in the way described 
in Sec. H imposes some restrictions on the kind of ob¬ 
servables that can be considered. Due to that, we con¬ 
sider here measurement operators that can be written as 
quadratic polynomials of the Majorana operators given 
in Eq. (3). 

Let us first consider the preparation of an arbitrary 
excited state of the XY Hamiltonian by applying U to 
the eigenstates of H[a]. For convenience, let us combine 
the internal parameters g and 6 of the XY Hamiltonian 
H, given by Eq. (8), in one variable ip = {g, S). According 
to Eq. (17) the eigenstates of the diagonal Hamiltonian 
H{(p)[a], are the computational basis states \k) while its 
corresponding energies are given by 

n—l 

Ek{(p) = Eo + '^e'l‘{ip), ( 68 ) 

1=0 

where ki are the binary components of k and the ener¬ 
gies €i are given in Eq. (18). Let us denote by I4'fe((p)) 
the eigenstate of E[{ip)[c\ corresponding to Ek{ip). As 
the product of matchgates U{ip) presented in Sec. HI 
maps H{ip)\c\ into the diagonal Hamiltonian E[{ip)[a\ [see 
Eq. (19)], we have 

\^k{^)) = U{^)\k). (69) 

Thermal states of the Hamiltonian E[[c\{ip) can also be 
generated by evolving the thermal state of the Hamilto¬ 
nian E[[a] with the unitary U{ip) [10]. Note that the ther¬ 
mal state of the diagonal Hamiltonian El[ip) [a] is given by 
p{ip,T)[a] = Jfc) (/cl, where Z{ip) de¬ 
notes the partition function Z{ip) = Efc Due 

to the fact that 

^ n— 1 

( 70 ) 

3=0 

the state p{ip,T)[a\ can be factorized as 

n—l 

p{ip,T)[a] = (^Pk{ip,T)[a]. (71) 
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The local density matrices pk{^,T) are diagonal and are 
given by 

T)[a] = Ukip, T) |0) (0| + T) |1) (1|, (72) 

where ak{p,T) = (1 + and bk{'^,T) = 

ak{(p,T). The thermal state of the Hamiltonian 
H{ip)[c] is obtained by evolving p((p,T)[a] under U{{p), 
that is 

p{v,T)[c] =U{ip)p{if,T)[a]U\ip). (73) 

Thus, thermal states and any eigenstates of the XY 
Hamiltonian can be prepared evolving the system that 
is initially prepared in a product state with the unitary 
U{ip). The final element of the circuit is a measurement 
of some observable on these states. As mentioned above, 
in order to be able to compress the matchgate circuit, the 
measurement operators are chosen as arbitrary operators 
A, acting on n qubits, which can be written in term of 
the Majorana operators Xj given in Eq. (3) as 

2n-l 

A = ^ Qj^kXjXk- (74) 

j^k=0 

In the following subsections we present the compressed 
quantum circuits running on m + 1 — log(n) + 1 qubits, 
corresponding to matchgate circuits whose outputs are 
the expectation values {A){ip,T) = ii {p{ip,T)[c\A} or 
{A){ip,k) = (4'fe((/9)|A|4'fc((,£))), where p{(p,T) (|4'fe((p))) 
denotes the Gibbs state (eigenstate) of the XY Hamil¬ 
tonian of n qubits respectively. We will then consider 
three examples of physically-relevant measurement oper¬ 
ators that can be written in the form given in Eq. (74), 
and will demonstrate how time evolution can be studied 
with an exponentially smaller system. 


A. Construction of the compressed circuit for 
thermal states 

Here, we consider the matchgate circuit that leads 
to the measurement of the observable A on the ther¬ 
mal state at temperature T. As shown above, the 
outcome of this circuit can be written as {A)(ip,T) = 
tr {U(if)p{(p,T)[a]U^(if)AY Using Eq. (5) and the ex¬ 
plicit expression of the observable A given in Eq. (74) we 
have 

2n-l 

{A){ip,T)= aj^k{j\R{‘p)S{(p,T)R^{(p)\k), 

j,k=0 

(75) 

where the 2n x 2n matrix S{ip,T) is defined as 
[S'((p, T)]^^^ = tr {—iXrXsp{p>,T)[a]}. Using the expres¬ 
sion for p{(p,T)[a] given in Eq. (72), one can see that 
the entries, [5'(v3, r)]r,s of S{(p,T) vanish unless (r, s) = 


{21, 21 + 1) or (r, s) = {21 + 1, 21), for some 1. Explicitly, 
we find that 

S{p, T) = Y tanh [(|2/) (2/ + 1| - |2; + 1) (2/|). 

m 

Let us define the positive operator p{Lp,T)[a\ = 
^ [I — iS{(p, T)] acting on m -I- 1 qubits. As we will see 
this operator corresponds to the initial state of the com¬ 
pressed circuit running on m -|- 1 qubits. One can easily 
verify that this state can also be written as 


^ n— 1 

p{‘P,T)[a] = ® l+y) i+y 

1=0 

bl{(p,T) \—y) { — y\^ , 


where |-|-y) and \—y) are the eigenstates of Y correspond¬ 
ing to the eigenvalues 1 and —1 respectively. Using 
Eq. (75) and the fact that R is orthogonal, the expec¬ 
tation value {A){ip,T) can now be written as 

{A){ip,T) = ntr [R{p)p{ip,T)[a]R^{ip)A] . (78) 

Here, we have defined the observable A acting on the 
Hilbert space of to -I- 1 qubits as 

2n-l 

A = iY <^j,k (lj) (^1 - \k) (j|) • (79) 

j,k=0 

With all that, the expectation value of any observable A 
[see Eq. (74)] can be measured using the following com¬ 
pressed circuit: 

(i) prepare the initial (to-I- l)-qubit state p{ip,T)[a], 
given in Eq. (77), 

(ii) evolve the state according to the orthogonal matrix 
R{(p) given in Eq. (52), 

(hi) measure the operator A, given in Eq. (79). 

Below, we consider three particular examples of phys¬ 
ically relevant quantities which can be measured with 
observables that are of the form given in Eq. (74) and 
therefore can be measured with a circuit as the one de¬ 
scribed above. Note that A can always be mapped to a 
local observable requiring at most 0{n) elementary gates. 


1. Magnetization 

We discussed in Sec. H that the XY Hamiltonian ex¬ 
hibits a phase transition as a function of the parame¬ 
ter g [cf. Eq. (8)]. As the magnetization shows an 
abrupt behavior at the critical value, pc, it can be uti¬ 
lized to observe this phase transition. In [7, 8] we de¬ 
rived compressed circuits to measure the magnetization 
of the Ising and XY model by simulating an adiabatic 
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evolution. With such a circuit, the magnetization of the 
groundstate for any value of g can be measured. Here, 
we derive a compressed circuit to measure the magneti¬ 
zation of both thermal states and arbitrary excited states 
of the XY Hamiltonian. In contrast to [7, 8 ] we will not 
simulate here the adiabatic evolution, but simulate the 
evolution of U(ip). 

The magnetization is the expectation value of an ob¬ 
servable M = which can be written as 

n ^-^3 3 

. n—1 

M = X2jX2j + l. (80) 

^ j=0 

Using this expression for M and Eq. (79), one can easily 
compute the operator M which needs to be measured in 
the compressed circuit. It is given by 

. n— 1 I 

M = --'£{ |2j) (2j + l|-|2j + 1) (2j|) = - (I 0 Y^), 

(81) 

where is the Pauli operator acting on the qubit. 
Using now Eq. (78) we have that 

(M)((^,r) = tr[i?((^)p((^,T)i?^(<^)(]l®U„)] . (82) 

Hence, the magnetization of the thermal state of the 
XY Hamiltonian, {M){ip,T), can be measured with the 
following compressed circuit running on log(n)-|-l qubits. 
First the initial m -I-1 qubit state p(<p, T) for the desired 
temperature [see Eq. (77)] is prepared. The gate R, de¬ 
composed as in Eq. (52), is applied and finally the last 
qubit is measured in y-direction. 

Eig. 1 shows the magnetization of the thermal states 
of the XY Hamiltonian as a function of g. The curves 
plotted correspond to two different system sizes, n = 8 
in Fig. 1 (a) and n = 128 in Fig. 1 (b). We have also plot¬ 
ted the magnetization corresponding to thermal states at 
temperatures between T — 0 and T = 0.9. In case T = 0, 
one observes that the magnetization exhibits a disconti¬ 
nuity at the critical point gc- This behaviour is due to a 
level crossing occurring between the groundstate and the 
first excited state, when Jordan-Wigner boundary con¬ 
ditions are considered (see Sec. H). This sudden change 
in the groundstate of the system, as a function of 5 , is 
reflected as a discontinuity of the magnetization of the 
system. Given that this behavior is due to the boundary 
conditions, it is a finite size effect and becomes negligi¬ 
ble as the system size grows. This can be seen when 
comparing figures (a) and (b). Note also that the tem¬ 
perature removes the discontinuity of the magnetization 
for all of the considered temperatures. This is expected 
as for any T > 0 and any value of 5 , the Gibbs state 
presents a mixture of the ground and excited states with 
some nonzero Boltzmann weights. Let us remark that 
this discontinuity is not due to the phase transition. In 
fact, if one considers only eigenstates with a given parity 
and momentum, the magnetization of the groundstate at 
T = 0 remains continuous at the critical point (see also 


[ 8 ] and references therein) and the phase transition can 
be witnessed as a discontinuity of the second derivative 
of the magnetization. 



FIG. 1: (Color online) The magnetization [see Eq. (82)] of 
the thermal states of the XY Hamiltonian as a function of g 
is shown. We have set 5 = 0.2 and considered system sizes 
of n = 8 qubits (a) and n = 128 qubits (b). Each curve 
corresponds to a different value of the temperature, which is 
chosen equally spaced between T = 0 (uppermost curve on 
the left side of the figures) and T = 0.9 (lowermost curve 
on the left side of the figures). At T = 0 the magnetization 
shows a discontinuity at the critical point (dashed line) as 
a consequence of the level crossing between the groundstate 
and the first-exited state. The discontinuity diminishes as 
the system size increases [compare figures (a) and (b)] and 
vanishes if T > 0. 



2. Quantum Quenching 


In [ 8 ] we derived a matchgate circuit and its corre¬ 
sponding compressed circuit to measure the number of 
kinks generated in a chain of Ising qubits, after the pa¬ 
rameter g is quenched from a large value, ^max, down to 
0. The number of kinks quantifies the number of spins 
that are not aligned in the X direction with respect to 
their neighboring spins. This number, normalized by the 
number of spins, can be quantihed hy i/ = 5(1 — (K)) [ 8 ], 
where the observable K is defined as 


K = 



n-2 

j=0 


n—2 

■^^^^^X2j+lX2j+2- 

i=o 


(83) 


We present now a circuit which can be utilized to mea¬ 
sure ^ as a function of the quenching time t, i.e., the time 
used to quench the parameter g, for thermal and excited 
states. As mentioned before, v corresponds to the num¬ 
ber of kinks only for the groundstate. However, as v can 
be also viewed as a measure of the correlations in the 
system, we investigate here its value for thermal as well 
as excited states. 

For convenience, we separate the circuit into two parts. 
In the first part the Gibbs state of the Hamiltonian 
^(5max)[c], p( 5 max,U)[c], is prepared. This is accom¬ 
plished by initially preparing the system in the state 
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p(,gmax,T)[a\ given by Eq. (71), and evolving it with 
the unitary [/(^max) presented in Sec. III. The corre¬ 
sponding compressed initial state and orthogonal gate are 
p( 5 max,T), given in Eq. (77) and i7(5max) [see Eq. (52)] 
respectively. 

The second part of the evolution corresponds to the 
quantum quenching. In this stage of the circuit, a prod¬ 
uct of matchgates UQ{t) simulates the quenching of the 
parameter g from (/max down to 0. We have provided the 
derivation of the unitary UQ{t), and its corresponding 
compressed gate Rgit) in [ 8 ]. Here, we simply recall the 
final decomposition of the matrix Rgit). The unitary 
UQ{t) consists of a discrete product of Trotter steps, ap¬ 
proximating the continuous change of the parameter g. 
Denoting the number of steps by L -I- 1, the error accu¬ 
mulated by this approximation scales as 0{LA'^), with 
The compressed gate Rgit) corresponding 
to the product of matchgates UQ{t), is given by 

L 

RQ{t)=l[Ro{t,l)Ri{t), (84) 

1=0 

where 

Ro{t, 1) = cos [2giA{t)] 1 + sin [2giA{t)] ho 
Ri{t) = cos [2A(t)] 1 + sin [2A(t)] hi, 

with gi = gmax (l - and 

n—1 

fT '0 = Yl (| 2 j -b 1 ) ( 2 j| - h.c.) 

n (86) 

/ri = ^(| 2 j + 2 ) ( 2 j + l|-h.c.). 

3 = 1 

Using Eq. (79), one can see that the observable K that 
has to be measured in the compressed circuit is given by 

n-2 

K = -^ ^ (|2j + 1) (2j + 2| - |2j + 2) (2j + 1|). 

(87) 

With all that, it is straight forward to write the expecta¬ 
tion value {K) (t) as the outcome of a compressed quan¬ 
tum circuit like the one given above [see Eq. (78)]. We 
obtain 

{K){t) = tr [i?Q(t)i7(5max)p(5max,T)[a]i?^(5max)77Q(0^] • 

( 88 ) 

Hence, the expectation value of the operator K can be 
measured as follows: 

(i) prepare the initial m + 1 -qubit state p((/max, T)[a], 
given in Eq. (77), 

(ii) evolve the state according to the orthogonal matrix 
i?Q(t)i?(gmax) given by Eq. (84) and Eq. (52), 

(hi) measure the operator K, given in Eq. (87). 


In Fig. 2 (a) and (b), we depict the correlations v, 
computed numerically using Eq. ( 88 ) as a function of g 
during the quenching for different total quenching times 
t. We have considered a value ^max = 10. Each curve 
depicted in the figures corresponds to a different total 
quenching time, chosen between tmin = 1 and tmax = 300, 
given in absolute units. The uppermost curve in both 
figures corresponds to a short quenching time, t = 1, 
while the bottommost continuous curve corresponds to 
the larger time t = 300. For comparison we have used 
the matrix R{g) to compute the correlations (for T = 0) 
that would occur in the system if the evolution were to 
be performed adiabatically. The curve computed in this 
way is plotted with a dashed line in Fig. 2 (a). 

The Kibble-Zureck mechanism establishes at T = 0 a 
relation iy{t) (x with p = 1 / 2 , between the num¬ 

ber of kinks and the quenching speed (1/t). We use the 
data plotted in Fig. 2 (a) to verify that such a mecha¬ 
nism is present in the Ising chain at T = 0 (see also [ 8 ]). 
In order to do so, we plot the logarithm of the corre¬ 
lations V at the end of the quenching evolution, i.e. at 
(/ = 0 as a function of the logarithm of the quenching 
speed. Due to the Kibble-Zureck mechanism, we expect 
a linear behavior of this function, with a slope p = 1/2. 
The data points have been fitted by a linear function in 
the region where a linear behaviour is observed, yield¬ 
ing p = 0.51 for a system size of n = 128 qubits [see 
Fig. 2 (c)], in agreement with the Kibble-Zureck mech¬ 
anism. The case where T > 0 lies out of the scope of 
the Kibble-Zureck mechanism and the relationship be¬ 
tween the correlations v (which at T = 0 correspond to 
the number of kinks) and the quenching speed is up to 
our knowledge unknown. In order to obtain some insight 
in the case where T > 0, we used Eq. ( 88 ) to compute 
numerically, for various temperatures, the correlations v 
present in a system initially prepared in a thermal state. 
As an example the values of In (ly) computed in the case 
where T = 2 are plotted in Fig. 2 (d) as a function of 
ln(l/t). In order to extract the proportionality factor p, 
we linearly fitted this function using the same interval of 
data points used for the case where T = 0 [see Fig. 2 (c)] 
and found that p = 0.03. In Fig. 3 we plot the value of 
p obtained in the same way when considering tempera¬ 
tures T G [0,2]. For low temperatures, there is a regime 
(T < 0.25) where the system still responds according to 
the Kibble-Zureck mechanism, as can be seen in Fig. 3. 
This behaviour can be understood as follows. For small 
enough temperatures only the groundstate and the first 
excited states have a non-negligible Boltzmann weight in 
the thermal state. Therefore, the evolution of the ther¬ 
mal state during the quenching evolution is determined 
by the evolution of the groundstate and the first exited 
states. The compressed circuit presented in this section 
can be easily adapted to simulate a quantum quenching in 
systems prepared initially in any eigenstate of the Hamil¬ 
tonian [26]. Doing so, we computed the ratio p between 
In(i^) and ln(l//:) for each of the first four excited states 
and found out that each of these values lie in the interval 
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[0.48,0.61] (for an example see Fig. 4 where we depict 
the logarithm of the correlations v present in a system 
initially prepared in the first and fourth excited states). 
Hence, to a good approximation, the dependence of the 
correlations v with respect to the quenching speed agrees 
with the Kibble-Zureck mechanism for low temperatures. 
However one can see that, as expected, the dependency 
of V on the quenching speed is drastically changed for 
higher temperatures (see Fig. 3). 



p 0.7 



0 0.5 1 0.75 


2 
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FIG. 3: (Color online.) The parameter p that corresponds to 
the exponent in the relation v oc is depicted as a fnnc- 
tion of the temperature T. Each point has been obtained 
by simulating a qnantum qnenching for a system of n = 64 
qubits. For each value of the temperature we have extracted 
p by fitting the correlations v measured in the system after 
the quantum quenching as a function of the quenching time t 
(see Fig. 2). Note that for low temperatures, the parameter 
p is close to the one for T = 0. This is expected as the corre¬ 
lations in each of the first four excited states scale similarly 
as the ones of the ground state. 


FIG. 2: (Color online.) The number of kinks, i/, as a function 
of the parameter g, occurring during the quantum quenching 
[see Eq. (88)] for a system size of n = 128 spins is depicted. 
The different curves correspond to various total quenching 
times ranging from t — 1 (uppermost curve) to t = 300 (con¬ 
tinuous bottommost curve). The number of Trotter steps was 
set to L = loot. The number of kinks for T = 0 and adiabatic 
evolution is depicted in (a) as a dashed line. Figures (c) and 
(d) show a standard linear fitting function of the logarithm of 
the correlations v at the end of the quenching, i.e. at gr = 0, 
vs the logarithm of the quenching speed (1/t). To fit the data 
points we have only considered the points depicted by filled 
circles, which correspond to long quenching times, although 
not long enough to correspond to an adiabatic evolution. The 
slope of the fitted function (represented by the dashed line) 
is the parameter p, corresponding to the exponent in the re¬ 
lation V oc We have obtained a value p — 0.51 for T = Q, 
and p — 0.03 for T = 2. 



FIG. 4: (Golor online.) The linear fit of the logarithm of the 
correlations u at the end of the quenching, i.e. at g = 0, vs 
the logarithm of the quenching speed (1/t) is shown. The 
simulated system size is of n = 128 qubits. The fitting was 
done in the same way as explained in Fig. 2. To obtain the 
correlations depicted in (a) and (b) we have considered a sys¬ 
tem initially prepared in the first and fourth excited state of 
the Ising Hamiltonian respectively. The corresponding slopes 
are p — 0.61 (a) and p — 0.48 (b). 


3. Correlations 

The compressed method of quantum computation con¬ 
sidered here can also be used to measure correlations 
(see also previous section). We focus here on those cor¬ 
relations which can be measured between two arbitrary 
qubits, say the and the fc*** . This correlation is quan- 
tihed by the expectation value of the observable 

Cj^k = + t • ■ ■ ^k—l^k — tt^2j + l‘^2k- (^^) 

As in the previous subsections, the observable which 
needs to be considered in the compressed circuit can eas¬ 


ily be computed. It is given by 

C,,fc = -i{\2j + 1) {2k\ - \2k) {2j + 1]). (90) 

Therefore, the correlations between the and the 
qubits of a thermal state is given by 

{Cj,k)i‘P,T) = ntr [R{(p)p{(p,T)[a]R^{ip)Cj^k], (91) 

where p{ip, T)[a], is given in Eq. (77) and R{ip) is given in 
Eq. (52). In Eig. 5, we depict the correlations occurring 
in a system of 64 qubits. In particular, the correlations 
between the central qubits, i.e. at position ^ and the 
qubit at position j is depicted as a function of j. To 
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compute the curves depicted in Fig. 5 (a), we have con¬ 
sidered a fixed value g = 0.8 and various temperatures 
between T = 0 and T — 0.9. The correlations decrease 
rapidly with distance and with temperature. Addition¬ 
ally, the correlations also decrease as the parameter g 
increases. This is shown in Fig. 5 (b), where each curve 
corresponds to a different value of g for T = 0. 



j 



j 


FIG. 5: (Color online.) The correlations of the form Cj^k 
[see Eq. (89)] between the jth qubit and the (^)th qubit as a 
function of j are depicted for n = 64 qubits (not all depicted). 
In (a) the dependency of the correlations on temperature is 
shown. Each curve corresponds to a different temperature 
whose values have been chosen to be equally spaced between 
T = 0 (corresponding to the uppermost curve) and T = 0.9 
(corresponding to the bottommost curve); the value of g has 
been fixed a.t g = 0.8. In (b) we show the dependency of the 
correlations on the value of (?, at T = 0. 


B. Construction of the compressed circuit for 
excited states 


In this subsection we extend the compressed circuit 
presented above to the case where the measurements are 
performed on arbitrary eigenstates of i/((p)[c], |'I'fe((p)) 
[see Eq. (69)]. In order to do so, let us define the operator 

n—1 

( 92 ) 

3=0 

such that \k) = Wk |0)®”. Using Eq. (69) we can write 
the eigenstates of H{ip)[c\ as |iE'fc((^)) = U{ip)Wk |0)®". 
The expectation value of the observable A measured on 
the state |'l>fe((/ 3 )) can then be written as 


(A)(<^, k) = tr \u{g3)Wk (|0) (0|)®” WkU\g3)A\ . (93) 


To obtain the compressed gate Rwk corresponding to 
the unitary Wk, we use a generalization of the procedure 
used to obtain the compressed gate Rbo in Eq. (60). That 
is, we write the conjugation of the Majorana operator Xk 
by the Pauli operator Xj as 

2n-l 

X,XkX,= ^ [Rx,]k,^i. (94) 

1=0 

Here, the matrix Rx^ is diagonal, with diagonal entries 
[Rxj]ii = I for Z < 2j and [Rxj]ii = — 1 otherwise. 


Thus, the compressed gate Rwk — n;: 
diagonal, with components 

1 /? 1 =/^ forr = I, 

I (—I)2fi(=i for r = 21 or r = 21 + 1. 

(95) 

The initial state of the compressed circuit can be com¬ 
puted in a similar way as p{g)[a\ in Eq. (77). We find 
pM = [cf- Sec. VA], with [^((p)]^^, = 

ti'{—iXrXs jo) (0|®"). This state can be written as [27] 

pH = ^ ® l+y)m (+2/lm • (96) 

Therefore, the expectation value of an observable A mea¬ 
sured on the excited state |'I'fe((/3)), can be written as 

{A){ip,k) = ntr [R{ip)RwkP[a]RwkR^i.^)^ ■ (97) 

The magnetization, correlations v, or correlations 
{Cij) of any excited state |4'/j((/)) can be measured by 
choosing as observables in Eq. (97) the operators M, K 
and Ci^j respectively. For example, the magnetization 
{M){ip,k) of an arbitrary excited state |4'/j((^)) is given 
by 

k) = tr [R{g,)RwkP[a]Ri;VkR^i^)Ym] ■ (98) 

As before, the compressed circuit can be directly read off 
this expression. 

In Fig. 6 the magnetization of each eigenstate of the 
XY Hamiltonian is shown. Each curve corresponding to a 
given eigenstate is obtained by computing k) as a 

function of 1/g for a fixed value of k and 5. Given that the 
number of curves grows exponentially with the system 
size n, we have considered here a small system size of n = 
8 qubits. Let us note that due to the level crossing at the 
critical point g^, the magnetizations have a discontinuity 
at this point, and the curves corresponding to different 
excited states cross each other (cf. Fig. 1). 


C. Time evolution 


The product of matchgates U recalled in Sec. HI can 
also be used to decompose the time-evolution operator 
V(t) = into a product of matchgates. In this way 

it is possible to simulate the measurement outcome of a 
time-evolved initial state with a matchgate circuit. In 
the following we briefly discuss the construction of this 
matchgate circuit and of the equivalent compressed cir¬ 
cuit. 

Let us first consider the initial state (|0) (Oj)®"; the 
output of the matchgate circuit described above is thus 
given by 


(A)(t)=tr [U(t)(| 0 ) {0\f-V\t)A 


(99) 


As in the previous subsections, we assume that the ob¬ 
servable A is of the form given by Eq. (74). Due to 
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efficiently generated with matchgates from the state 


(| 0 ) ( 0 |) 


0n 


To obtain such a circuit it is sufficient to in¬ 


FIG. 6: (Color online.) The magnetization for different eigen¬ 
states of the XY Hamiltonian as a function of 1/g and for 
S = 0.2 and n = 8 is shown. Due to the boundary condi¬ 
tions, several magnetization curves exhibit a discontinuity at 
the critical point (see Sec. II), as indicated by a vertical line 
in the figure at the value = (1 -I- d)”^ = 0.83. 


Eq. (19) we have V{t) = UW{t)W, where W(t) = 
g-zff[a]t diagonal Hamiltonian H[a] as a function of 
Pauli operators is given in Eq. (70). It follows that W{t) 
can be decomposed as the following product of match- 
gates 


Wit) = 


n 

i=o 




( 100 ) 


In order to obtain the compressed gate associated 
to W{t), we write H[a] as a quadratic polynomial 
in the Majorana operators using the Jordan-Wigner 
representation given in Eq. (3) and obtain H[a\ = 
^ {x 2 jX 2 j+i — X 2 j+iX 2 j)- From this expression, 
and following the procedure described in Sec. II one can 
see that the compressed matrix Rw{t) corresponding to 
W (t) is given by 


n—1 

Rw{t) = ^ ( _ 
1=0 


cos{ejt) sin(ejt) 
sm{€jt) cos{ejt) 


( 101 ) 


we 


In the same way as in the previous subsections, 
rewrite the expectation value {A){t) as 

{A)it)=tT{[RRwit)R^Pi^RR^{t)R^A), ( 102 ) 

where Pin = <8 l+n)™ (+t/L, and the observable A can 

be computed from A using Eq. (79). Thus, we have that 
a compressed quantum circuit comprising of the prepa¬ 
ration of the initial state pin, the evolution of the system 
under the action of the orthogonal matrix RRw{t)R^ 
and the measurement of the observable A can be used to 
compute the expectation value {A){t) given in Eq. (99). 

The matchgate circuit considered above can be eas¬ 
ily generalized to include any initial state that can be 


sert at the beginning of the matchgate circuit described 
above the matchgates that prepare the desired initial 
state. The compressed circuit has to be modified ac¬ 
cordingly by the set of compressed gates corresponding 
to the additional matchgates. 

A more general set of initial states can be consid¬ 
ered in the case where the restriction to a logarithmi- 
cal space classical computer is removed. In this case, 
one can consider an arbitrary initial state of the form 
Pin = Po ® • 0 Pn-i, where each pi is a density ma¬ 

trix of a single qubit. The initial state that one has to 
consider in the equivalent compressed circuit is given by 
Pin oc [al — iS{p)], where S{p) is defined by [S'(p)]j = 

tr {—ixjXkp), where the Majorana operators Xj are given 
in Eq. (3). The real parameter a has to be chosen in 
such a way that pin is positive-semidefinite. Note that 
the preparation of pin requires the classical storage of all 
the components of S'(p), which is of dimension v?. 

Let us note that in [8] we presented a compressed quan¬ 
tum circuit that can be used to simulate the time evolu¬ 
tion of a XY chain of qubits. In contrast to the circuit 
presented here, in [8] the compressed quantum circuit was 
constructed starting from a Trotter decomposition of the 
time-evolution operator leading to certain errors in the 
estimation of {A) it) that depend on the simulated time 
t. Here, by using the matrix U that exactly diagonalizes 
the XY Hamiltonian, {A) (t) can be computed exactly for 
any value of t. 


VI. CONCLUSION 


We derived compressed quantum circuits to measure 
physically interesting properties of the ID-XY-model. In 
order to do so, we used the decomposition into a product 
of matchgates of the unitary U that diagonalizes the XY 
Hamiltonian, and the fact that any eigenstate or thermal 
state of the system can be prepared by evolving a prod¬ 
uct state according to U [10]. Due to the equivalence be¬ 
tween matchgate circuits and compressed quantum com¬ 
putation [6], the corresponding matchgate circuits can 
be compressed as long as the measurement operator is 
quadratic in the Majorana operators. That is, there is 
a compressed quantum circuit that simulates the match- 
gate circuit, running on exponentially less qubits. In the 
cases considered here, we have additionally shown that 
the compressed circuit requires not only less qubits, but 
also less gates than the matchgate circuit that it simu¬ 
lates. Let us finally mention that the results presented 
here can be extended to any quadratic Hamiltonian. The 
reason for that is that the matrix that diagonalizes any 
of these Hamiltonians can always be decomposed into a 
polynomial number of matchgates [12]. 
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Appendices 

A. REORDERING TRANSFORMATIONS 

The aim of this appendix is to provide a more de¬ 
tailed description of the product of matchgates corre¬ 
sponding to the various reordering operations mentioned 
in Sec. IIIC [see Eq. (48)]. Let us recall that each 
of these transformations can be decomposed into ele¬ 
mentary transformations which map the pair of 

fermionic operators {aj,ak) into the pair {ak,aj). The 
transformations Sj^k has the associated matchgate Sj^k 
given in Eq. (50). Therefore, the product of matchgates 
corresponding to an arbitrary reordering transformation 
S can be easily obtained from the decomposition of S in 
terms of transformations Sj^k- 

A convenient way of providing such a decomposition 
is by associating to the transformation S an operation 
that permutes the components of the vector fj. intro¬ 
duced in Eq. (21) in Sec. III. Eor example, let us 
consider the transformation Sj^k that maps the pair of 
operators {aj,ak) into the pair {ak,aj). This trans¬ 
formation has an associated unitary matrix Sj^k that 
transforms an arbitrary state |4'[a,/x]) into a new state 
Sj^k |^[a, m]) = [see Eq. (21)], where = /Zfe, 

n'k = fij and /i[ = /i;, ^ j,k. Therefore, we can 

associate to the transformation Sj^k the operation that 
permutes the components j and k of the vector fi. In the 
subsequent sections we apply this procedure to decom¬ 
pose the transformations for 0 < s < m and iSeog 
in terms of elementary transformations of the form Sj^k- 


Reordering transformation 

The reordering transformation 5®°® is the first step of 
the transformation U given in Eq. (48). Therefore, the 
unitary matrix 5'®°® associated with this transformation 
acts directly on the initial state of the quantum circuit. 
The quantum initial state can be associated to fermionic 
operators Uj according to Eq. (21). This association es¬ 
tablishes a particular ordering of the operators aj that is 
described by a vector = (0 ,..., n — 1). 

The reordering transformation 5®°® precedes the Bo- 
goliuvov transformation B, given in Sec. Ill A. The lat¬ 
ter maps the operators {ai}i=o,...,n-i into the operators 
{&i}i=o,...,n-i) by mixing the pair of operators aj with 
a-j, for —j = n — j, and ao with a^. In order that 
these operators act on consecutive qubits [according to 
Eq. (21)], the matrix 5'®°® has to transform the vector 


into some vector A®°® whose components are paired 
in the form (j, — j) and (0, |). This condition still al¬ 
lows for multiple choices of the vector A^°®. In this work 
we have chosen a particular ordering with the aim of 
obtaining a simple expression for the compressed circuit 
associated to the matrix 5'®°®. This ordering is given by 
the vector 

= (O, 2, -2,... 21, -21, ..., ^ - 2, -^ + 2, 

--l,-- + I,...,2/-l,-2/ + I,...,l,-lj. (103) 

It is simple to verify that A^°® can be obtained from 
A*-°^ by permuting successively the components ctq = f 
and aj, where 

77- 77- 

aj = - + (-iyj, for j = 0,...,--I. (104) 

Given that the permutation of two components, e.g. (Jq 
and CTj is associated with the fermionic swap gate -S'o-o.o-j, 
it follows that the unitary matrix corresponding to the 
transformation 5®°® is given by the following product of 
matchgates 

■2 

n (105) 

i=i 


Reordering transformation for 1 < s < m — 1 

The reordering transformation acts after the 

transformation and before the transformation 

presented in Sec. IIIB. Recall that the transfor¬ 
mation maps a set of operators into 

the set by mixing pairs of operators 

whose indices differ solely on the (m — s — 1)*^ com¬ 
ponent. The matrix associated with the transfor¬ 
mation acts on neighboring qubits in case where 
the association between qubits and fermionic operators 
is given by the vector A*-®^ introduced in Sec. IIIB. Let 
us recall that the components of this vector are given 
by = n^s and = n^s + 2""“®“\ where n^s is 

given in Eq. (42). Therefore, the reordering transforma¬ 
tion 5i®i has to be associated with the operation that 
maps the vector A*-®”^^ into the vector A*-^^ 

In order to see which components of A^®~^^ have to be 
permuted to obtain the vector Ai®\ we now write down 
the binary decomposition of the elements of these two 
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vectors. Using Eq. (42) we have that 


(5-1)' 

^21 

— [^m—2 • 

• ■ s+1 

0 ^m—s 

^m — s—1 ■ 

■ ■ ^0 

(s-l)- 
^2/ + l 

= [^m—2 • 

• ■ s+1 

1 ^m—s 

^m — s—1 ■ 

■ ■ ^0 

\(s) 

= [^m-2 • 

• ■ s+1 

^m — s 0 

^m — s—1 ■ 

■ ■ ^0 

1 

= [^m—2 • ■ 

• ■ s+1 

^m—s 1 

^m — s—1 ■ 

■ ■ ^0 


(106) 

One can see that only half of the components have to 
be reordered and therefore, j permutations are required. 
Let us consider, for example, the 2/*^'' component. For a 
value of I where Im-s = 0, we have that 
thus this component does not have to be reordered. The 
same happens with the (21+1)**' component when Im-s = 
1. The only components that have to be permuted are 
the {2iy^ and the (21' + I)*** , in case where I is such 
that Im-s = 1 and I' = I® 2™“^. Note that this is where 
the binary decomposition of I' is equal to the one of I for 
any bit except the (m — s)*** . This can be written in a 
compact way as follows. Two components zi(r)('*^ and 
Z 2 (r)^'^^ of the vector must to be permuted if they 


take the form 




U(r)(®)' 

— ■ • 

■ 5 s 51 ? '^m— s — 1 5 ■ • 

■ -+0,0], 


[z2(r)(®^j 

— 35 ■ ■ 

• s— 15 ■ ‘ 

■ ■ Uo, 1], 


(107) 

for some 0 < r < q — 1 with a binary decomposition 
given by [r] = , ro]. 

In the following we define the set 


flgCs) = |(^Zi(r)(®\2:2(T’)^'*^)} 


r=0....f-l 


(108) 


where zi(r)^®^ and Z 2 (r)^'''^ are given in Eq. (107). Note 
that the set fisc.) groups the pair of indices that have to 
be permuted with each other in order to transform the 
vector into A*'®^ Using this definition, it follows 

that the product of matchgates associated to the trans¬ 
formation 5*-®^ is given by the product of fermionic swap 
gates 


= n ( 109 ) 


Reordering transformation 

The reordering transformation 5^**^ acts after the 
transformation iSgog and before the transformation 
[see Eq. (48)]. In the same way as in the previous sub¬ 
sections we associate with the transformation 5^**^ an 
operation that permutes the components of the vector 
/j,, whose components indicate which qubit is associated 
with which fermionic operator. 


Given that the reordering transformation iSgQg inverts 
the reordering produced by iSsog, after this transforma¬ 
tion, the association between qubits and fermionic op¬ 
erators is given by the vector A*-**^. However, in order 
for the transformation to act on pairs of operators 
with consecutive indices, it is necessary that the associa¬ 
tion between qubits and fermionic operators is given by a 
vector A[see Sec. Ill Bj. Therefore, the transformation 
^(0) niust correspond to the transformation of the vec¬ 
tor A*-**^ into the vector A^*\ Using Eq. (42) we can write 
the binary decomposition of the components of these two 
vectors as 


\(o) 

'^21 

v(0) 

^ 2+1 


( 1 ) 


^ 2+1 


[^m — 2 ■ ■ • 

0 

0 

[^m — 2 ■ ■ • 

^0 1] 

[0 l-iji—2 • 

• ■ ^ 0 ] 

[1 l-fji—2 • 

• ■ ^ 0 ] 


( 110 ) 


In order to obtain the permutations required to transform 
A^°^ into A‘'^\ we can use the results of the previous sec¬ 
tion. Comparing the Eq. (110) to Eq. (106), we see that 
the transformation can be obtained as a concatena¬ 
tion of the transformations 5^®\ with s from m — 1 to 
1. It follows that the unitary matrix associated with the 
transformation 5^**^ is given by the product 

s<o) = ( 111 ) 

where the gates for 1 < s < m — 1 are given in 
Eq. (109). 


Reordering transformation 


The reordering transformation acts after the 

transformation jr(iTi-i) IIIBj. Therefore, the 

unitary matrix corresponding to this transforma¬ 

tion yields the output state of the circuit. In the follow¬ 
ing we see how to obtain the permuting operation of the 
vector /X associated with 

Let us stress that in both the input and the output 
states, the association between qubits and fermionic op¬ 
erators Uj and Cj respectively, must be given by the same 
vector A*'*’^ Note however that the matrix yields 

states where the association between qubits and the 
fermionic operators is given by the vector 

[see Sec. IIIBj. That means, using the notation intro¬ 
duced in Eq. (21), that the output state is a state of the 
form |■i/;[a;^™“*\ A*^™^]). However, the operators 
are related to the operators Cj according to Eq. (41), 
i.e., = cj', with [j] = [jm_i,..., jo] and [f] = 


[jO) • ■ ■ ) jm—l] ■ 


Therefore, a state 



equal to the state 


vI'[c,a'] 


where the binary decompo- 
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sition of the components of the vector A is given by 



bo) ■ ■ ■ ) jm—l]■ 


( 112 ) 


From the discussion above, it follows that the ma¬ 
trix must be associated with the operation that 

transforms the vector A into the vector A^*^^ with 
= [jm-i) • ■ •) Jo]- This transformation can be 
achieved as follows. Let us consider a transformation 
that maps the vector A^ to another vector A in such a 
way that the and {m — t — 1 )**' binary components 


of the elements of 


are permuted. That is 


[jO? • ■ • ? Ji —15 jm—t — 1 ijt+li * • ■ : jm—t — 2 ^ jm—t: • ■ • 7 Jm—l] 
Applying these transformation for 0 < t < tmaxi with 
tmax = L^/ 2 J to the vector A' yields the vector A*'®^ 
Denoting by for 0 < t < tmax the product of 
matchgates associated with each of the aforementioned 
transformations, it follows that the matrix can be 
written as 


_ 2 ^(imax) . . . y( 0 )_ (113) 


We now describe how to obtain the decomposition into 
matchgates of each of the unitaries . From the discus¬ 

sion above we have that is associated with a trans¬ 
formation of the vector A^ that permutes two components 
(r) and (r) if and only if they are of the form 



yi\r) 

= [r 771-3, 

■ ■ • , 0, '^771—3—t, ■ • ■ 

,n, 1 , 



n-i,- 

••Uo], 



y 2 ^ (r) 

= [^771-3, 

■ ■ • ,'^771 —A —t, ^,'^771 — 3 — 1, ■ • ■ 




n-i,- 

■•Uo], 

(114) 

for any 0 < 

VI 

Defining the set 



HtU) = 1 [j/i ^ 

(r),V^(r)]| , 

J r=0....5-l 

(115) 

where (r) and y^'^ 

(r) are given in Eq. (114), 

it follows 

that the matchgate can be written as 



= n (116) 

U,k)en^(t) 


Hamiltonian. Moreover, we provided a decomposition 
of these matrices as a product of matchgates Sj^k- In 
this section, we provide the compressed orthogonal ma¬ 
trices that correspond to these products of matchgates. 
The compressed gate corresponding to Sj^k, is given by 
[see Eq. (67)]. Therefore, a direct way to obtain 
the compressed gates corresponding to the product of 
matchgates A®”®, for 1 < s < m — 1 and S’!™) 

is to use the same decompositions given in Eq. (105), 
Eq. (Ill), Eq. (109) and Eq. (113) respectively, in terms 
of the matrices Rsj,k^ instead of the matchgates Sj^k- 
One can easily verify that these decompositions require 
O [poly(n)] matchgates. However, due to the symmetries 
of the problem, it is possible to rewrite the compressed 
gates corresponding to entire reordering steps with at 
most O [log(n)^] elementary gates, as we show within 
this appendix. 


Compression of the reordering transformation 

for 1 < s < m — 1 


Here, we consider the construction of the compressed 
matrix Rg(s ), corresponding to the product of matchgates 
that has been given as a decomposition in term of 
fermionic swap gates in Eq. (109), for 1 < s < m — 1. 
As we pointed out above, we can decompose Rsu) in 
term of matrices Rs k using the decomposition given in 
Eq. (109), that is 


Rsi‘)= n ^s„k- ( 117 ) 

(j,fc)eng(^) 


Using the explicit form of the matrix Rs^ ^ given 
in Eq. (67), and the fact that none of the indices 
j occurs more than once as an element of a pair in 
Eq. (117) can be written explicitly as Rg{s} = 

1 + E(i.fc)en(o) ( \j) (fc| + \k) (j| - |j) (j| - \k) {k\ 


im- This expression can be written in a more compact 
way by using the binary decomposition of the indices 
(j, fc) in the set i.e., using Eq. (107) leads to 


With this transformation, we conclude the construc¬ 
tion of the product of matchgates corresponding to each 
of the reordering transformations. In the following ap¬ 
pendix, we provide the compressed orthogonal gates cor¬ 
responding to each of these product of matchgates. 


B. COMPRESSION OF THE REORDERING 
TRANSFORMATIONS 


r-l 


Rsisy= l+^|r)(r|®(|01)(10| + |10)(01| 


r =0 


-jOl) (01|-|10)(10| 

= (jOO) (00| + |11)(11| 

+ | 10 ) ( 01 | + | 01 )( 10 | 


) s— Ij-m— 1 


^ Im 


) s— l,m— 1 


0 1 , 


(118) 


In Appendix A we presented unitary matrices acting 
on quantum states that are associated to the reorder¬ 
ing transformations used in the diagonalization of the 


where the identity operator acts on all qubits except on 
the (s — 1)*'' and the (m — 1)*'' . Note that Rgi^) is the 
swap gate between the qubits s — 1 and m — 1 , i.e. a 
single two-qubit gate. 
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Compression of the reordering transformation 

Due to Eq. (Ill) we have that the compressed gate 
i? 5 (o) can be written as 

i?5(o) = i?5(i) • • •-R5(m-i), (119) 

where the matrices Rgis), for 1 < s < m— 1 are two-qubit 
gates, given in Eq. (118). Therefore, Rsm is a product 
of m — 1 elementary gates. 

Compression of the reordering transformation 5^"*^ 


Generalizing this equation to the case where we have f — 1 
factors we obtain 


Rs...= E 10(^1 +El C’i) (Ci+ll + lo’f-l) (cqI 1 Glim 

i =0 J 

(123) 

where we have defined E = {o'j}j=o,...n-i- Using the bi¬ 
nary decomposition of the indices aj, it is possible to 
rewrite the previous expression in the following more 
compact way 


To construct we need the compressed gates 

Rrp(t}, corresponding to the product of matchgates given 
in Eq. (116). Given the similarity between the definitions 
of and the matrix Rrp(t} can be constructed 

analogously to Rs(s)- We obtain 


RgBog — 


(|00) (00|o,m-l + |ll)(ll|o.m-l))' 


+ |10) (01|o.m-l ® ^ + I01) (10|o.m-l ® BA 


(124) 


Rtm =(| 00 ) (001 + 111 ) (111 

+ |10) (01| + |01)(10|)^ _^_^G1. 


( 120 ) 


where the operators A and B act on qubits 1,..., m — 2 
and are given by 


Hence, R^i*) is the swap gate between the qubits t and 
m — 1 — t. Due to Eq. (113) it follows that 


Rg(m} — i?7’(t„ax) ■ ■ ■ Rt( 0 ) ■ (121) 

As each matrix Rx(t) is a two-qubit gate, the implemen¬ 
tation of Rgim.) requires tmax = [m/2j two-qubit gates. 


k 

^max 

-4 = E 1 ^^ (fcmax - k\ , 

k^O 

k —1 

'^max 

B = E! 1^ + 1) (^1 + |9) (fcmaxi ) 
k =0 


(125) 


Compression of the reordering transformation 

From the decomposition of S'®”® given in Eq. (105), it 
follows that its corresponding compressed matrix can be 
written as Rgsag = n/=i BSfTQj ? where the indices (7j 
are defined in Eq. (104). Using the explicit expression of 
the indices aj and the matrices Rg^ ^ given in Eq. (67), 
it is possible to express Rgsag in a more compact way. 
To do so, note that the product of two matrices Rao,aj 
and Rao.a-k takes the form 

= (1 - ko) (ctoI - \(Jj) {(7j\- IcTfe) (cTfcl 
+ ko) (Cjl + \aj) (CTfcl + |CTfc) (ctoI ) G Im- (122) 


with /cmax = n/A — 1. Note that A = 
is a product of to — 1 single-qubit gates. The 
operator B acts on the computational basis state 
as B \j) = \j 01) where © is the addition mod¬ 
ulo 2"^“^ and can be decomposed as a product of 
TO — 2 controlled operations. More explicitly, we have 
that B = [A(—2) (X„_ 3 )] ... [A(™-2).-.2(Xi)], 

where {Ok) denotes the single-qubit controlled op¬ 

eration O acting on the qubit fc, and the ii,... are the 
control qubits. Given that any controlled operation act¬ 
ing on TO qubits can be decomposed into a product of 
0{m) elementary gates [19], it follows that Rgs^g given 
in Eq. (124) can be decomposed into a product of 0{w?) 
elementary gates. 
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/ \ \2l,2l+l 

Using Eq. (43) to write the 


> as 

21 , 21+1 

x<=) x(“) x<=) x(^) 

^ 2 i ^21 ''21 + 1 ^21 + 1 


transformed operators < x 


do+l) .^(n + 1) 


x( = ) 


x<=) 

^21+1 
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of the operators 


C) AA 


<»)’ 

‘21 


X 


(s) 

21 + 1 
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